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Abstract 

The newly observed D* sJ family containing 1)^(2317), £> sJ (2460) and D sJ (2632) 
attracts great interests. Determining their structure may be important tasks for both 
theorists and experimentalists. In this work we use the heavy quark effective theory 
(HQET) and a non-relativistic model to evaluate the production rate of D*j(2317) 
from the decays of ■0(4415), and we find that it is sizable and may be observed at BES 
III and CLEO, if it is a p-wave excited state of £> s (1968). Unfortunately, the other 
two members of the family cannot be observed through decays of charmonia, because 
of the constraints from the final state phase space. 

PACS numbers: 13.25.Gv, 12.39.Hg, 12.38.Lg, 14.40.Lb 

1 Introduction 

The recently observed exotic mesons D*j(2317), D s j(2460) and -D s j(2632) [Q seem to 
constitute a new family of mesons which are composed of charm and strange flavors. The 
mesons possess spin-parity structures of + , 1 + and + respectively. This new discovery 
draws great interests of both theorists and experimentalists of high energy physics. Some 
authors 2,. suppose that D*j(2317) and D s j(2460) are the chiral partners of the regular 
D s and D*, while D s j(2632) may be a radially excited state of _D*j(2317). They may also 
be considered to be p-wave excited states of D s |5]. Alternatively, many authors suggest 
that they can possiblly be four-quark states or molecular states [U E] • The most peculiar 
phenomenon is that in some experiments the three resonances are observed with clear 
signals pQ, whereas not by other prestigious experimental groups. One would ask if the 
observed resonances actually exist or the background was misidentified as a signal. It is 
noted that similar situations exist for pentaquarks [Sj. The goal of the research is to help 
designing experiments which can help clarifying the mist. 

The key point is to experimentally explore the resonances and find a convincing expla- 
nation why they are observed in certain experiments, but not in others. However, before it, 



one needs to design certain experiments to confirm the existence of D*j(2317), D s j(2460) 
and -D s j (2632) and determine their hadronic structures. As aforementioned, there are 
several different postulates. Measurements may tell which one is more realistic. 

Because of the constraint of final-state phase space, observing a final state which 
involves any of the exotic states can only be realized via decays of higher excited states in 
the ifj family. From the data-booklet [7j, we can see that the lowest excited state which can 
offer sufficient energy to produce D* sJ (23\7) + Z) s (1968) is -0(4415), but still not enough 
for D*j(23l7) + 5^(2317). However, since D* sJ (2317) is a 0+ meson and £> s (1968) is 
a CP meson, a careful analysis on the total angular momentum and parity indicates the 
decay of -0(4415) — * D*j(2317) + Z) s (1968) is forbidden. Moreover, if only considering 
the central value of ■0(4415), the phase space is not enough for D*j(2317) + Z) s (1968) + tt 
and -D*j(2317) + Z)*(2112) which could be produced via pure strong interaction and the 
only possible mode is the radiative decay -0(4415) — > D*j(2317) + Z) s (1968) + 7. That is 
a decay with a three-body final state and is an electromagnetic process where a p-wave 
is necessary to conserve the total angular momentum and parity. This observation tells 
us that the corresponding branching ratio must be very suppressed and is a rare decay. 
Recently, Barnes et al. [H] suggest to observe Z)*j (2317) via the process 

0>(4415) -► D* J (2317) + J D*(2112). 

The advantage is that it is a strong decay with a two-body final state, therefore the 
amplitude may be large, but meanwhile, by the central values, m +ro n »™„i > 4415 

MeV, thus this reaction can only occur via the threshold effect and would suffer from a 
corresponding suppression. If it is of a larger rate (we will estimate it later in the work), 
the decays 0(4415) -» D* i7 (2317) + 5 s (1968)+_7r and 0;(4415_) -> D* sJ {23l7)+D_ s (1968) +7 
can also be realized via secondary decays of Z)*(2112) — ► £> s (1968) + tt and Z)*(2112) — ► 
l) s (1968) + 7 and these are the dominant modes over the direct decay modes 0(4415) — > 
D* sJ (2317) + D S (1968) +vr and 0(4415) -> D*j(2317) + D s (1968) + 7 which are three-body 
decays. 

The picture is that the charmonium -0(4415) dissolves into a cc pair and both c and 
c are free and on mass shell, and the soft gluons emitted from cc can excite the physical 
vacuum to create a pair of ss. The process of ss pair creation is quantitatively described 
by the quark-pair-creation model (QPC)|H1 E3- Then the ss join the corresponding c 
and c to compose charmed mesons. Indeed, the creation process is fully governed by the 
non-perturbative QCD effects, thus the rate is not reliably calculable so far and can only 
be estimated in terms of models. In this work, we use QPC model |HJ El to evaluate the 
rates of ^(4415) -> £>*j(2317) + D*(2112) and the direct decay 0(4415) -» D* J (2317) + 
-D s (1968) + 7 where a photon is emitted during the process. 

In this work, we consider the transitions of -0(4415) — ► D* sJ (23l7) + Z)*(2112) and the 
subsequent observable modes -0(4415) -> D*j(2317)+-D*(2112) D*j(231t)+D s (1968) + 
7 (D*j(2317) + 5 S (1968) + it). We also calculate the ratio of the direct radiative decay 
process -0(4415) — > D*j(2317) + l) s (1968) + 7 which is not produced via the resonance 
DJ(2112). 

The key point is how to evaluate the hadronic matrix elements. Here we must adopt 
suitable models to do the job. 
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Z)*j (2317) and Ds all are heavy mesons, therefore one can expect that the heavy 
quark effective theory (HQET) applies for evaluating the hadronic matrix elements. For 
a completeness, we keep the l/m c corrections in the formulation, however, it is obvious 
that such corrections are practically negligible in the concerned case, so that we do not 
really include them in the numerical calculations. As a check, we employ a non-relativistic 
model to re-evaluate the hadronic matrix elements and compare the results obtained in 
the two approaches. 

To obtain the concerned parameters and testify the applicability of the model, we 
calculate the branching ratios of 0^(4040) — ► 13 M + £)(*) an d D s + D s . By fitting data, 
we determine the vacuum production rate of the quark-pairs in HQET. Moreover, when 
using the non-relativistic model, we also need to determine the concerned parameters in the 
wavefunctions. More concretely, there are several decay channels in -0(4040) with c and c 
in the final states (_D°(*).D (*), .D^MjD^M), and their branching ratios are experimentally 
measured. Actually, in HQET, the only free parameter is the rate of quark-pair creation 
from vacuum, i.e. 7 g , then one mode is enough to fix it. We can check the obtained 
model and the parameter by applying them to evaluate other modes which have also been 
experimentally measured. Our numerical results respect the pattern determined by the 
experiments. For -0(4415) more channels are available, that is D S D S (or D+D~), etc. 
We may naively consider that the production of D S D S in 0(4415) — > D s + D s is somehow 
related to -0(4040) — * £)W+jDW, then all the parameters obtained from decays of -0(4040) 
can be applied to study decays of 0(4415) while assuming the parameters are not very 
sensitive to the energy scale. 

In both HQET and non-relativistic model, we derive the formulation for the branching 
ratio of ^(4415) -> D*j(2317) + D*(2112) and obtain the final numerical results. We 
also formulate the direct process -0(4415) — > D*j(2317) + Z) s (1986) + 7 which is the only 
channel allowed by the phase space if neglecting the threshold effects. Even though these 
results with the aforementioned approximations cannot be very accurate, one expects that 
the order of magnitude of the calculated result would be right. 

If the exotic states _D*j(2317) is of the 4-quark structure as suggested 013, in the 
production process at least three pairs of quarks are created from vacuum, and the final 
state would involve more quarks and anti-quarks, thus the integration over the final-state 
phase space would greatly suppress the rate. By our rough numerical evaluation, at least 4 
orders suppression would be resulted for the decays, if the exotic meson D*j(2317) is a four- 
quark state. Thus by measuring the branching ratio of -0(4415) — > D*j(2317) +D* (2112), 
we may judge (1) if the exotic meson £>*j(2317) indeed exists, (2) what quark structure 
it possesses. 

This work is organized as follows, after the introduction, in Sect. 2, we formulate 
the decay rates of ^(4415) -► D*j(2317) + Z)*(2112) and direct process ^(4415) -> 
D*j(2317) + Z) s (1968) + 7. In Sect. 3, we present our numerical results along with all 
the input parameters. Finally, Sect. 4 is devoted to discussion and conclusion. Some 
detailed expressions are collected in Appendix. 
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2 Formulation 



The QPC model about the process that a pair of quarks with quantum number J = 
++ is created from vacuum was first proposed by MicujU] in 1969. In the 1970s, QPC 
model was developed by Yaouanc et al. ^01 EEH EE2 EE] and applied to study hadron 
decays extensively. Recently there are some works EH to study QPC model and its 
applications [lfij . In the QPC model, the interaction which represents the mechanism of 
a pair of quarks created from vacuum can be written as |15j 

Svac = g iq J d^Xtpglpg, With C = g iq i)q^q, (1) 

where g Iq = 2m q r y q . ~f q is a dimensionless constant which denotes the strength of quark 
pair creation from vacuum, and can only be obtained by fitting data. m q (q = u, d, s) 
are the masses of light quarks. In the non-relativistic approximation {T2|, the interaction 
Hamiltonian Q can be expressed as 

U vac -> nZc = E / d Pq dp q [3 lq 5 3 (p q + P g) E(i, i; m, -m|0, 0) 

i,j m 

xyr(Pa-Pg)(xI m ^ ) l M(Po,s)d+(p q ,S > ), (2) 

where i and j are SU(3)-color indices of the created quarks and anti-quarks; s and s 1 are 
spin polarizations; ipo = (uu + dd + ss)/y/3 and (u>o)ij = Sij for flavor and color singlets 
respectively; xT 1S a triplet state of spin, y™ is a solid harmonic polynomial corresponding 
to the p-wave quark pair. 

2.1 The transition amplitude of -0(4415) in QPC model. 

In this work, we study the strong decay ■0(4415) —* Z)*(2112) + D*j(2317) and the direct 
radiative decay V(4415) -> 7 + 5 S (1968) + £>*j(2317). With the QPC model, during 
these transitions charm-quark and antiquark from ■0(4415) combine with the ss created 
from vacuum to form final state particles. The Feynman diagrams of these transitions are 
depicted in Fig. ^ 

The transition matrix element of -0(4415) -> D*(2112) + D* sJ (2317) is 

T stron g = (d*(2112)D* sJ (2317) \H vac (x) |0(4415)). (3) 

For the direct radiative decay -0(4415) -> 7 + l) s (1968) + D*j(2317), the transition 
matrix element reads as 

T dir = (5,(1968)^(2317)7^ J dxdy[C vac (x)C em (y)MU15)), (4) 
where C em (y) is the electromagnetic interaction Hamiltonian and have the following form 

£em(y) = ±y / d 4 x* 7 ^A»(y). (5) 
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Figure 1: (a) is the Feynman diagram of strong decay ^(4415) — > D*(2112) + D*j(2317). 
(b)-(f) are the Feynman diagrams for direct ^(4415) -> 7 + S s (1968) + D*j(2317). 

where the sign ± corresponds to charges of c and c respectively. Considering the weak 
binding approximation, |?/>(4415)) can be expressed as 

|^(4415)) -» NV(0)cfc\0), (6) 

where ^(0) is the wave function at origin and e denotes the polarization vector. N is the 
normalization constant. 

It is also noted that for decays ^(4040) — ► D^D^*\ the Feynman diagrams are similar 
to that in Fig.l a. 

2.2 Evaluation of the hadronic matrix elements in HQET. 

(i) The strong decay ^(4415) -> 5*(2112) + D*j(2317). 

The diagram in Fig. ^ (a) involves the q-meson-Q vertices. In refs. |17l I18j . the 
effective Lagrangian for these vertices has been constructed based on the heavy quark 
symmetry and chiral symmetry 

C HL = h v {w-d)h v -[gx{H + S)h v + H.c]+g'[Tr{HH) + Tr{SS)\, (7) 

where the first term is the kinetic term of heavy quarks with fh v = h v ; H is the super-field 
corresponding to the doublet (0~, l - ) of negative parity and has an explicit matrix repre- 
sentation: H = ~2'{P^ ti — -P75); P and P*^ are the annihilation operators of pseudoscalar 
and vector mesons which are normalized as (0|P|M(0 - )) = \/Mh, (0|P* m |M(1 - )) = 
s[M^€»\ S is the super-fields related to (0+,l+) and S = ^[^7^75 - P ]l X = & 
(q = u,d,s is the light quark field and £ = e t ' , here we only take the leading order as 
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The effective Lagrangian in Eq. contains only the leading order for the coupling of 
meson with quarks. We may also include the l/mg corrections in the expressions. The 
heavy-light quark interaction Lagrangian given in ^Hj is 



2g 2 



A- 



Q{i$ — m Q )Q — j^Qlu—Q^—^, 



(8) 



where Q = (b, c) and ip = (u, d, s), m Q is the heavy quark mass. We can obtain the l/m c 
corrections from two aspects. The first comes from the quark wavefunction ^5] 



Q(x) 



-tm_ v-x 



1 + 



2m, 



h v (x), 



- v^v ■ D, 



(9) 



the l/m Q correction is obtained by replacing h v with (1 + H~/>^/2m Q )h v in (J7J). Secondly, 
the superfields H and S in (J7| should also receive l/m g correction. Falk et al. [2U| 
presented the changes as 



H 



2m, 



2m, 



Then, we can include the l/m Q corrections in l[7|). 



Now, let us write down the transition amplitude for the decay of ■0(4415) — ► -D*(2112) + 



2^(2317) as 



M(i/>(tel5) -► D*(2U2) + £>* 7 (2317)) 



*(0) 



Tr 



x(M A +i> A )U 



c 



i> c -i> A l 2 - r "~ 



-9u 



-g^M, 



2*(0)g 2 g Ia VM B M c 



3^M 



i A 



M A mi 



-M 3 A + M A (p A - Pc ] 



m. 



M A M 2 C 



B 



[(p A /2-p B ) 2 -m 2 s ][(p A /2-p c ) 2 -ml}' 



(10) 



where p A , p B and p c represent the four momenta of ^(4415), D*j(2317) and Z?*(2112); 
M A and e A are the mass and the polarization vector of ■0(4415); Mb and M<j are the 
masses of two produced mesons; g is the coupling constant of Q-meson-q vertex which is 
given in literature JlZj- It is noted that by the central values 



m. 



™ oJ(2112 ) > 4415 MeV ' 



thus the process can only occur through the threshold effect. The resonance ■0(4415) 
has a total width F A . Considering the distribution, we adopt the typical Gaussian form 
suggested by the data group [7], and set the lower and upper bound for the integration 
of final phase space as M A — 5 < M < M A + 5 and the delta-function guarantees the 
energy-momentum conservation. 
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Finally we obtain the width 

r(0(4415) -» ££(2112) + D*j(2317)) 

< M - M "> 2 W (id 



where 



2(r A /2)^ 



IPj = m » (12) 



E B = VMl + pl, E C = ^MI + V %- S = lM-f, f3 = W%U\. (13) 

For the indirect subsequent decays ^(4415) -> ££,(2317) + ££(2112) -» ££,(2317) + 
££(1968) + 7 and -0(4415) -► ££,(2317) + ££(2112) -» ££,(2317) + ££(1968) + vr, the 
rates are obtained as 

r nd (V>(4415) -> ££,(2317) + ££(1968) +7) 
= r(0(4415) -> ££,(2317) + ££(2112)) x BR(D*(2112) -> ££(1968) + 7), (14) 

r nd (^(4415) -> ££,(2317) + 5,(1968) + vr) 
= r(0(4415) -> ££,(2317) + ££(2112)) x 5i?(££(2112) -> ££(1968) + vr). (15) 

(ii) The direct radiative decay 0(4415) -» 7 + ££(1968) + ££,(2317). 

This process is much more complicated than the strong decay depicted in (i). 

By the QPC model a pair of ss quarks is created from vacuum and the underlying 
mechanism is the soft gluon exchanges which excite the vacuum sea. The momentum of 
the light quark pair created from vacuum is small and the photon hardly has possibility 
to be produced from light quark. Thus we can ignore the contribution of Fig. ^ (d) and 
(f) for the direct 0(4415) -> 7 + ££(1968) + ££,(2317) process. 

The amplitude for radiative decay 0(4415) -> 7 + £> s (1968) + ££,(2317) includes 
several pieces. For Fig. E(b), 

M {b) = Q C ^J j^^^ iS2 ( q )[v( P 2,s 2 )0^u( Pl ,s 1 )i (16) 

{b) = -ys9y/Mc-t , \ n ~ ■ 9is ~L , L 1 — ^—g^M^ 



i>j2 — ft — m 
for Fig. ^ (c) , the transition amplitude reads as 

T f dq 



ic-i>A^-^s la ~i> B + j> A /2-ft-m s 

fa- 




M ic) = -QJ—J j^A 1s MMp*,s 2 )OMu(pi,si)), (17) 
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o (c) = h 



# ~i> A /2 - m, 
i 

' i> A l 2 -i> B - m > 



Aa/ 2 ~^b ~ m - 



■9u 



9\[M> 



B 



where pi, pi and s\, are the momenta and spin projections of the charm quark and 
anti-charm quark, and the following relations hold 

Pl+P2=P A = (M A ,0), Pl -pa = 2q = (0,2q), £ J dq|^ SlS2 (q)| 2 = 1. 



*1,S2 



Using the method of |21j . we have 




(2vr)3/2 



i>s 1S2 (q)«(P2, s 2 )O w w(pi, s x ) 



^ ; r _ )T!3/2 -^ lS2 (q)rr[(M A 0« + {C#\fl + + M M • 0«)1±^(-^)](18) 

where of = O (i) | q =o and d«) = g|rO W |q = o, and i denotes (b) or (c). 

Dissociation of the charmonium into cc can be well described by the non-relativistic 
model where the wavefunction at origin ^(0) corresponds to the binding effect. (|16l) and 
((T7|) can be further expressed as 



M 



(b) 



Tr 



■9u 



B 



i> A j2 — ft — m 



f -t A /2-m B * ls -j> B +j> A l2 



m. 



M 



(c) 



X9i 



Tr 



ft —j) A /2 - m, 



i A l 2 -i> B - m > 



gVMB(M A +p' A )tA 



(19) 



(20) 



where p A , p B , p c and k correspond to the four momenta of ■0(4415), D*j(2317), D s (1968) 
and photon respectively; e& is the polarization vector of the emitted photon. 



The decay width for 0(4415) — ► 7 + D s (1968) + D*j(2317) radiative decay is expressed 



as 



1 



n 



( d 3 p t 



T y2ir) 4 5 4 (M A -p B - Pc -k)\M [b) +M [C) \ 



6M A J Y \ (2Tr) 3 2Ei 
In next section, we carry out the multiple integration to obtain numerical results. 



(21) 
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2.3 Evaluation of the hadronic matrix elements in the non-relativistic 
model. 



As discussed above, for a comparison 1 we are going to employ a non-relativistic model 
i.e. the harmonic oscillator model to repeat the calculations made in terms of HQET. 
Application of such model should be reasonable in this case. 

(i) Strong decay 0>(4415) D*(2112) + £^(2317). 

We calculate the ^(4415) -> D* (2112) + D*j(2317) by using QPC model in the non- 
relativistic approximation. The decay width is 

r(0(4415) -► 5*(2112) + D*j(2317)) 

r ™ n _E B Eo\k\ , a _r ( M ~ M -r dM ^ m 



(1 - 0)y/2ir T a Jm a -8 



/ 2rr — 2JMa exp 



2(r A /2) 2 



and the concrete expression of J2i s \Mi s \ is collected in Appendix. The definitions of 8 and 
(3 in the expression are exactly the same as those given in ea.(|13|). 

(ii) For direct radiative decay 0(4415) -> 7 + 5,(1968) + D* sJ (2317). 

Following the traditional method [13], the matrix element of radiative decay ■0(4415) — ► 
7 + Z) s (1968) + D*j(2317) in nonrelativistic approximation can be written as 



non 
vac 



(VdM, 4; P4, 4)^(2317) (P'i, P3. 4)* 7 ( k > e(k))|T[W 

•Wem]|^V(4415)(Pl> s l) P2> S2)) 

=t^e x^uxltt x^x? / n n 

in J o=l 6=1 

X<5 3 (pi + P2 " P^)<5 3 (p' 2 + P4 " P D > 3 (P3 + P4)^ 3 (p'l + P 3 " P DsJ ) 

x (1, 1; -n Dsj , n D JO, 0) (1, 1; n, -n|0, 0)^ n (p 3 - p 4 ) 

1 1 . , , 1 . 1 

X^iPl ~ 2 P ^/mP2 - -pJ^ DsJ (Pi - ^P DsJ ,P 3 - 2 Pd sJ> 

X^ S (P2- 2Pd s 'P4- 2 P ^)(°l 6 Pi d P^Pi 6 Pi a k y rf^yiT^c^tx) 

x&pA^pJ )' ( 23 ) 

where T and C correspond to the flavor and color factors in this transition; x's are the spin 
wave functions; p^, p Ds and p D are three-momenta of ■0(4415), Z) s (1968) and Z)*j(2317); 
(ftp, (p D r and tp Ds are the harmonic oscillator wave functions of 0(4415), D*j(2317) and 
-D s (1968) respectively. 

1 There is another reason to employ the non-relativistic harmonic oscillator model. If D*j(2317) is of 
four-quark structure, the HQET no longer applies and the only model we can use for the multi-constituents 
structure is the harmonic oscillator model. Therefore a comparison of the results in the model with that 
obtained by HQET is indeed meaningful. Namely, HQET is believed to be applicable in this case, thus 
consistence of the results obtained in the two approaches can confirm applicability of the harmonic oscillator 
model. Then we can use it to calculate the production rate if Z)*j(2317) is of four-quark structure. 
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(b) 



Neglecting some technical details, finally one can derive the decay amplitude of Fig. ^ 



M {b 0(U15) -> 7 + D* sJ (2317) + A,(1968)) 



9 7rV4 

1 



let 1 \ 2/3 1 



: (2vr) 4 



(ip 2 exp -- J Rl(2p 2 r--i?c-(p c + 2p fl +2p 2 



-i?|(2p 2 + Pl 



xy 1 nB (-2p 2 -p B )y 1 n (-2 Pfl -2p A +2p 2 )w(p 2 ,s 2 ) 7 ^(p c + Pfl -p A +p 2 ,s 2 /)e^(k) 
-105 + 210i$p= - 84 ^ + 8^p^ 



12V35 



where 5 is a spin factor, ^(0) is the wave function of ■0(4415) at origin. The indices A, B 
and C are for 0(4415), D* sJ (2317) and £> s (1968) respectively. 

With the same treatment, we also obtain the amplitude M( c ) °f Fig. ^ ( c ) ; an d for 
saving space we keep its expression in Appendix. 



2.4 A rough estimation of the production rate of D*j(2317) in ^(4415) 
decays, if D*j(2317) is of a four-quark structure. 

There have been some works which suggest that -D*j(2317) is of a four-quark structured, 
the situation would be completely different. We draw a possible Feynman diagram 
in Fig. 2, and one can notice that three quark-pairs are created from vacuum. As more 
particles are produced, the final state phase space would greatly reduce the rate. 



V>(4415) 




Figure 2: The Feynman diagram describing the production of -D*j(2317) in 0(4415) in- 
clusive decay considering the four quark structure of D*j(2317) in the QPC model, where 
ellipsis denotes the diagrams for other possible quark combinations. 

The inclusive transition matrix element can be written as 

<c, S , g ,g 1 2?; j (2317)|r[7C?(^lW(x 2 )^(x3)]|^(4415)). (25) 

Thus the decay width reads as 

r(V>(4415) -» £>*j(2317) + c + s + q + q) 
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x |M(0>(4415) -> Djj(2317) + c + s + g + g)| 2 . (26) 

Two points are noted: First this amplitude cannot be evaluated in the framework of 
HQET, but only in the non-relativistic model because of the complicated quark-structure; 
Secondly Fig. 2 depicts an inclusive process -0(4415) — > _D*j(2317)+4 free quarks, consider- 
ing hadronization, observable processes can only be 0>(4415) — > D*j(2317) + -D S (1968) + tt 
and -0(4415) — > D*j(23l7) + Z) S (1968) + 7 where qq annihilate into a photon. As discussed 
above, such direct processes are much suppressed. 

Because the inclusive decay ^(4415) — > D*j(2317) is related to multi-body final states, 
the calculation is very complicated. The multi-integration over the phase space is very 
difficult, even in terms of the Monte-Carlo method. Generally the rate is proportional to 



a ~ (7 9 ) 3 



which is a remarkable suppression factor. 



47T 

J2nf 



Thus if -D*j(2317) is of a four-quark structure, one can expect that the inclusive decay 
width of -0(4415) — > Z?*j(2317) is at least four orders smaller than the corresponding value 
if _D*j(2317) is a p-wave excited state of the regular D s meson. 



3 Numerical results 

(l)Determination of the concerned parameters in the two approaches, 
(a) The parameters for the transition in HQET. 

The coefficient G = yg 2 , which is introduced in the transition amplitude, is obtained 
by fitting the decay width of -0(4040) — * DD. In Appendix, we present the formula for the 
decay rate of charmonium into two charmed pseudoscalar mesons in HQET. For calculating 
G, the value of ^(4040) (0) is obtained by fitting the experimental data of 0^(4040) — ► e + e~ 
HU- We get ¥^(4040) (0) = 0.101 GeV 3 / 2 and G = 12.3 GeV" 1 . In ref.[T7|, the value of 
g 2 is obtained as g 2 = 4.17 GeV -1 , thus we obtain y q = 2.95. Yaouanc et al. used 
to employ the harmonic oscillator to evaluate such processes, and they got 7^ 3 [Tl?] . 
which is very close to the value we obtain with HQET. However, it is also noted that y q is 
purely a phenomenological parameter and its value may vary within a reasonable range, 
for example, in their later work, Yaouanc et al. took y q to be 4 instead, when they fitted 
datafTS], 

Since there are not enough data to determine 7 S , we adopt the relation 

7s = lq/VS, 

for later numerical computations. 
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(b) The parameters in the non-relativistic model. 

In this scenario, the non-relativistic approximation is taken and the expression is no 
longer Lorentz invariant, the relevant parameters may be somehow different from the 
values in HQET, especially the value of 7^ which corresponds to the vacuum creation of 
a quark pair. However, as pointed above the values obtained in these two approaches are 
very close, thus we can use 7^ = 2.95 for later calculations. 

Using the experimental results of V>(4040) DD, 0(4040) -> D*D* and ^(4040) -> 
DfD~ decays |23] 2 , we obtain all the relevant parameters which are needed for later 
numerical computations, For the readers' convenience the relevant formulations |12j are 
collected in Appendix. With all the information, we obtain the values of R's in the har- 
monic oscillator wave functions as: R\ = 6.00 GeV~ 2 , R 2 D = 5.25 ± 0.22 GeV" 2 , R 2 D , = 
6.70 ± 0.67 GeV~ 2 and R 2 Dg = 5.20 ± 0.58 GeV" 2 . 

However, R D * (2317) an d Rd* corresponding to D* sJ (2317) and D* cannot be obtained 
by fitting data, because such decay modes do not exist for the sake of phase space of final 
states. As priori assumed, D* sJ (2317) is a p-wave excited state of D s [H], therefore the 
difference in R of Z)*j(2317) and D* is due to the L ■ S coupling which is proportional to 
l/m c . Furthermore, due to the heavy quark symmetry, the difference in R of D*(2112) and 
D s is also of order l/m c , thus R 2 D , ( 2 317) ~ ^D* ~ R"b s are employed in the calculations of 
V>(4415) -> D*(2112) + D* sJ (2317) and direct decay 0(4415) -► 7 + D s (1968) + £>*j(2317). 
It is believed that this approximation is reasonable for estimating the order of magnitude 
of these transitions. 

(2) Our numerical results for the D*j (2317) production. 

In the calculations of 0(4415) -> D* (2112) +D*j (2317) and the direct decay ^(4415) -> 
7 + D s (1968) + Z)*(2317) by the two approaches, we employ the following parameters as 
inputs: M^ 415) = 4.415 GeV, M D ± = 1.968 GeV, M D * = 2.112 GeV, M D * j(2317) = 
2.317 GeV [Zj. By fitting the data of -0(4415) — > e + e~ which is available at present, we 
obtain ^(4415) (0) = 0.088 GeV 3/2 . 

We now present the numerical results obtained with the two approaches in Table. 2. 





I 


II 


Br(0(4415) -► D* (2112) + D* j (2317)) 


9.16% 


9.58% 


Br (0(4415) -»■ 7 + £> s (1968) + D*j(2317))(ind) 


8.63% 


9.03% 


Br(i;(U15) -► it + D s (1968) + D* J (2317))(ind) 


5.31 x 10" 3 


5.56 x 10" 3 


Br(0(4415) -► 7 + L> s (1968) + L>* J (2317))(dir) 


8.46 x 10" 5 


2.29 x 10" 5 



Table 1: Columns I and II correspond to the numerical results obtained in HQET and the 
non-relativistic model respectively. 

2 The BES measurements of the inclusive charm cross section at 4.03 GeV are |23|: a D a + ago = 
19.9 ± 0.6 ± 2.3 nb, a D + + a B - = 6.5 ± 0.2 ± 0.8 nb and a D + + a D - = 0.81 ± 0.16 ± 0.27 nb. Considering 
the relation [T3]: r(D*°D*°) : T{D*°D° + D°D*°) : T{D°D°) » 1 : 7 : 9, we obtain the following- 
decay widths of V(4040): r(^(4040) -» DD) = 2.97 ± 0.68 MeV, T(V'(4040) -» D* D*) = 26.73 ± 
6.13 MeV and r(V>(4040) -> D+D~) = 1.55 ± 0.69 MeV. 
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4 Discussion and conclusion 



It is obvious that the newly discovered D s j family may be very significant for better 
understanding of the hadronic structure and low energy QCD. The members of the family, 
D*j(2317), D s j(2460), D s j(2632), all have positive parity, so that they cannot fit in an 
s-wave cs(cs) structure. The literatures suggest that they may be p-wave excited states, 
namely chiral partners of D s , D* etc. or four-quark states as well as molecular states. It 
is necessary to look for a more plausible way to determine their configurations, i.e. design 
an experiment(s) to clarify the picture. At least we would like to find an experiment to 
judge (1) if such states indeed exist, (2) their quark configuration (p-wave excited states 
or four-quark states). 

To have a larger production rate, it is reasonable to look for .D*j(2317) via strong 
decays of higher excited states of charmonia. The most possibly available charmonium 
is 0(4415). Since -D*j (2317) mesons have positive parity, the decay mode of ■0(4415) —* 
D*j(2317) + Z) s (1968) is forbidden and if considering the central values of the masses of 
the concerned particles and constraints from phase space of final states, only ^(4415) — > 
D*j (2317) + D S (1968) + 7 is allowed. This direct radiative decay must be much suppressed 
as discussed in the introduction. 

Barnes et al. jS] suggested to observe decay 0(4415) ->■ D*j(2317) + L>*(2112) which 
can occur via the threshold effects. The consequent decays 0(4415) — ► D*j(2317) + 
L>*(2112) -> D* J (2317)+5 s (1968)+7andV(4415) -> L>*j(2317)+Z)*(2112) D*j(2317)+ 
ZX,(1968) + 7r can be observed. Even though such processes may only occur via threshold 
effects and should be suppressed, it is noted that m n , + rn is only slightly 

above 4415 MeV, one can expect that the suppression is not very strong. 

In this work, we carefully study the production of D*j (2317) in the decays of ^(4415) 
and evaluate its production rate. The processes are realized as the charmonium ?/>(4415) 
dissolves into a cc pair which then combines with ss created from vacuum due to the 
non-perturbative QCD effects and constitute two mesons. The first step is determined by 
the wavefunction of -0(4415) at origin and the light-quark-pair creation is described by 
the QPC model For evaluating the hadronic transition matrix elements, we employ 
two approaches, i.e. HQET and the non-relativistic model. Our final numerical results 
achieved in the two approaches confirm this allegation as they are reasonably consistent 
with each other. 

To guarantee the plausibility of the results, we obtain all necessary parameters by fit- 
ting data. However, it is understood that there must be some errors from both theoretical 
and experimental aspects, and the parameters should have some uncertainties, especially 
the vacuum creation rate of the light quark pair. Thus the real rates may be within a 
range around the values estimated with the input parameters and theoretical approaches, 
the order of magnitude should be correct and trustworthy. 

For a comparison, we have also evaluated the transition rate of the direct radiative 
decay 0(4415) — ► D*j(2317) + D s (1968) + 7 in the same approaches and find that the 
resultant rate is two orders smaller than that through the intermediate state Z)*(2112), 
even though it is realized via the threshold effects. 
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We find that even though the threshold effects suppress the production rate of ^(4415) — 
D*j(23l7) + jD*(2112), it is still sizable if D*j(2317) is a p-wave excited state. If so, it 
can be observed in the future experiments of BES III, CLEO and maybe at Babar or 
even LHC-b. However, as our calculations indicate that if D*j(2317) is of a four-quark 
structure, its production rate is much more suppressed and cannot be observed in decays 
of charmonia. 

Unfortunately, in such decays, one can only expect to observe D*j(2317), but not 
the two other members of the new family. However, once the existence and structure 
of -D*j(2317) are definitely confirmed, we have reasons to believe existence of the other 
two. Moreover, we can have more knowledge on the hadronic structure and may design 
experiments to testify the other two. We are looking forward to new experimental results 
to clarify this theoretical problem. In a recent work, some authors calculate the 
decay rates of D*j(2317) and D s j(2460). They claim that their results prefer the ordinary 
cs(cs) quark-structure for the mesons. However, a decisive conclusion must be drawn from 
a deterministic experiment (s), and D*j(2317) + D*(2U2) suggested by Barnes et al. as 
well as subsequent observable modes D* J (2317) + 5 a (1968)+7, D*j(2317) + 5 a (1968)+7r 
would provide an ideal possibility to make this judgement. 
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Appendix 

(a) In the Eq. the J2i s \Mu\ is 



„3/2 „3/2 „5/2 

tx A n B n c 
147456 v / 35vr?? 11 /2 



cxp 



k R A (R B + R\ 



i?4 



.3 1,6 



48a & (2C + l)rfk 



+8a V (2C(C + l)??k 2 + 3) k 6 - 336a 3 (2( + l)7? 2 k 4 + 12a V (l4((C + l)v^ 2 
+8V3 + 3) k 4 - 420a(2C + l)r?k 2 + 42a 2 r?(lOC(C + l)??k 2 + 8\/3 + 3^k 2 
+105(2C(C + l)^/k 2 + 9)] - 84r]R\ [ - 16a 3 (2£ + l)?7 2 k 4 + 4a 4 r/ 2 (2C(C + l)??k 2 + 3)k 
-40a(2C + l)r?k 2 - 4a 2 ??(2C(C + l)r/k 2 - 4^3 + 2l)k 2 + 15^2C(C + l)r/k 2 + 7 
+1680r/ 2 i? 2 1 [6C(C + l)f?k 2 -4a(2C + l)r/k 2 + 2a 2 r/(2C(2C + l)r/k 2 +3)k 2 + 15 

-6720r ? 3 (2C(C + l>/k 2 +3)], (27) 



where 



R 2 a 



r2 a + r2 b + r2 c 



r2 a + r2 b + R c ' 



i + C- 
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(b) The concrete expression of M. (c) is 



JW W (»>(4415) - 7 + »:(2317) + D,(1968)) 

i ,s m l^)l&f''^\-^(J-) 2 ' 3 1 fry* 

V^357V 7T J 9 vr 1 / 4 L 3 V2vr/ v / 2£fc 1 
x | d Pl exp [ - i J Ri(2 Pl ) 2 - ±R%(p B + 2 Pc + 2 Pl ) 2 - \r 2 c {2 Vi + p c ) 2 ^ 

x^r" B (p s +2 Pc +2 Pl )^i n (-2 Pc -2p 1 )fi(p B + Pc + Pl ,5 2 ) 7 XPi,S2')^(k) 
-105 + 210g|p] - 84gjpj + 8g|pj 
12v/35 



(c) The amplitude of "0(4040) decay into two pseudoscalar mesons in HQET is 

2ig. g 2 m a JM A M B M c #(0) 
_M(V(4040)-P + P) = A I ^2 2 12 ^-(Pb-Pc)- (29) 



3 [(PA/ 2 -Ps) 2 - m . 



3 J 



(d) In ref. ^2j , the authors gave a general expression for calculating decays of ^(4040) 
DD, 0(4040) -» + £>£>* and ^(4040) -> L>*L>* 



T(^(4040)) = Ck 3 N 2 (k 2 ), (30) 

where C means a spin-SU(3) factor corresponding to the particular channel under con- 
sideration (C = 1/3 for DD, C = 4/3 for DD* + DD* and C = 7/3 for D*D*), k is 
the three-momentum of the final particles in the CM frame of 0(4040), and N 2 (k 2 ) is a 
normalization factor and has the following expression 

^ (fc2) = 4^0^ [ ^ /2(4e)eXp( - e)]2 ' (31) 
where L^J 2 is a Laguerre polynomial, and £ = k 2 R 2 /6. 
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